Vibronic mechanism of high- Tc superconductivity 
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The dispersion of the in-plane Cu-0 bond-stretching LO phonon mode in the high- Tc super- 
conducting cuprates shows strong softening with doping near the zone boundary. We suggest that 
it can be described with a negative electronic dielectric function that results in overscreening of 
inter-site Coulomb interaction due to phonon-induced charge transfer and vibronic electron-phonon 
resonance. We propose that such a strong electron-phonon coupling of specific modes can form a 
^ ■ basis for the phonon mechanism of high-temperature superconductivity. With the Eliashberg theory 

using the experimentally determined electron dispersion and dielectric function, we demonstrate the 
possibility of superconductivity with the order parameter of the d k 2_ k 2 symmetry and the transition 
temperature well in excess of 100K. 

PACS numbers: 74.20.Mn, 74.20.-z, 74.72.-h, 74.20.Rp 
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I. INTRODUCTION 



The mechanism of the high-T c superconductivity (HTSC) in the cuprates |gj remains elusive, in spite of extensive 
experimental and theoretical efforts. The majority view on the mechanism is to consider magnetic interactions as the 
main driving force However, the HTSC recently observed in MgB2 [gj, graphite-sulfur composites Q, and the 
, n-type infinite layer cuprate Sro.gLao.iCu02 || cannot be explained by magnetic mechanisms, since these compounds 
t-H ' have little spin fluctuations in strong contrast to the cuprate superconductors. There is no reason to reject, then, the 
possibility that the HTSC in the cuprates also shares a similar non-magnetic mechanism. In this paper we discuss 
a phonon mechanism of HTSC based upon the overscreening of the inter-site Coulomb interaction, in light of recent 
experimental results on the cuprates. In this mechanism, unlike the BCS theory, the relevant phonons are not the 
long wave acoustic phonons but the zone-edge optical phonons that induce inter-site charge transfer. 

The HTSC cuprates are doped Mott insulators. Even though the antiferromagnetism disappears with only about 
2 % of hole doping, strong antiferromagnetic spin fluctuations are observed by neutron scattering, nuclear magnetic 
resonance, and other methods. For this reason spin fluctuations have been considered to be the principal mechanism 
of the HTSC However, the intensity of spin fluctuations does not correlate, or even anticorrelate, with T c 

O ■ as the hole concentration is changed. For instance the spin fluctuations measured by the nuclear relaxation of Cu in 
Tl2Ba2Cu06-,5 are the same for the sample with T c = 85 K and the overdoped sample with T c — 0. Even when 
the composition is the same, the spin fluctuations of L^-zBa^CuC^ films prepared by epitaxial growth that have T c 
of 47 K were found to be smaller than those of the bulk with T c of 30 K Jiy] . 

In this paper we propose a phonon mechanism based upon the anomalous screening of the inter-site Coulomb 
interaction in a highly correlated electron system leading to strong pairing. This mechanism does not compete 
against the magnetic mechanism, and could achieve HTSC alone or with a magnetic mechanism through synergetic 
effect. There is a large volume of literature that show coupling of superconductivity to the lattice and phonon 



Jll| , |l2 13 FJ. In particular the in-plane Cu — O bond-stretching FO phonon mode was observed by neutron inelastic 
scattering to show strong softening with doping near the zone boundary along the Cu-0 bond direction [ ll| Jl5| , ^6[ . 



This mode induces charge transfer between Cu and O, and thus couples strongly to the charge [|T^,^7 18 . Recent 
neutron scattering measurements on YBa2Cu30y_5 (YBCO) suggest that the frequencies of the FO phonons are 
strongly softened with doping near the zone-boundary, while the TO phonons are not |l9|j2C|]. We propose that this 
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occurs due to the phonon-induced charge transfer that results in formation of a vibronic state |2lJ and a negative 
electronic dielectric function as will be discussed in section II. 

Doped holes in the cuprate superconductors are highly correlated due to the on-site Coulomb interaction on Cu ions 
and resultant spin fluctuations. Holes are strongly renormalized due to these interactions, and as a result many physical 
properties of the cuprate oxides in the normal state show anomalous behavior. For example, very unusual temperature 
dependence has been observed for the Hall J2^] and Seebeck coefficients J2^] and the pseude-gap excitations [g4| have 
been observed above T c . The quasi-particle band structure determined by a ngular resolved photoemission is very 
different from those calculated with the LDA due to the correlation effect |2q| . 

Because of these strong correlation effects at present it is very difficult to calculate the electronic response function 
accurately from the first-principles. Instead, in the present paper we take a phenomenological approach based upon the 
information obtained from the experimental data. We start with the model electronic band structure that agrees with 
the angle-resolved photoemission experiments p5fl , and estimate the effective interaction using the phonon dispersion 
determined by the inelastic neutron scattering experiments p9| , p0| , p6[ . We then calculate the superconducting order 
parameter and the superconducting transition temperature on the basis of Eliashberg's formulation p7[] . We show 
that the order parameter is of the (4-2 _ fc 2 -wave symmetry p8| and the superconducting transition temperature can 
exceed 200K. Because of the d-symmetry the present phonon-mediated pairing mechanism does not compete against 
the spin-fluctuation mechanism. It is possible that both mechanisms operate in the cuprate, and the weight depends 
upon the doping level. The issue of the in-plane anisotropy in the phonon dispersion []19| , po|j2q ] and the effect of 
the possible local spin/charge stripe fluctuation will not be included in this paper, and will be discussed in future 
publications. 

In section II we explain the idea of overscreening in the cuprate superconductors. In section III we derive the 
equation with the charge kernel for the calculation of the order parameter and T c . In section IV we determine 
the energy band renormalized by the correlation effect by utilizing the experimental results of the angle-resolved 
photoemission and the dielectric function using the results of the neutron scattering measurements. In section V 
we determine the kernel of the equation by using the model obtained in Section IV. We then numerically solve the 
equation and obtain the superconducting order parameter and estimate the value of the superconducting transition 
temperature. The section VI is devoted to discussions. 



II. OVERSCREENING MECHANISM AND THE EFFECTIVE INTERACTION BETWEEN CARRIERS 

The starting point of this research is to recognize that the electron-phonon (e — p) interaction in such a strongly 
correlated electron system can be very different from that in the conventional metals. In particular, since the carrier 
density in the cuprate superconductors is relatively low and the charge dynamics is strongly correlated with the spin 
dynamics, the medium-range Coulomb interaction is not fully screened. For instance the optical reflectivity does not 
saturate as in the standard Drude model even below the nominal plasma frequency |29[] . Thus the system has strong 
dielectric interactions, unlike in the usual metals. In particular, since phonons modulate the covalent bonds they 
induce charge transfer between ions, and therefore local polarization. In this paper we suggest that this effect can be 
described with a negative electronic dielectric function. 

Sometime ago Tachiki and Takahashi proposed an overscreening mechanism of phonon-mediated superconductivity 
pO| , pT|j32| . To explain this mechanism let us consider a static case of the overscreening effect. We write a staggered 
external electric field as D(q) and the induced staggered electric charge polarization as P e z(q), and define the electronic 
charge susceptibility \ei (q, 0) by 

47rP e *(q)=Xe*(q,0)£>(q) (1) 
Then, combining an electromagnetic relation D(q) = E(q) + 47rP e ;(q) = e e ;(q, 0)E(q) with Eq.(l), we have 

' 1-Xei(q,fj). (2) 



e ei (q,o) 

According to the Kramers-Kronig relation, we have 



,oo j 

X e ;(q,0)=2/ du)~p el (q,uj), (3) 
Jo w 



where p e i (q, uj) is the spectral intensity of the electronic charge fluctuation given by 



p e ;(q,cj) = ~Im[ , J (4) 
7T e e i{(\,uj) 
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The spectral intensity is always positive and therefore we have an inequality from Eqs. (2.3,4), 



- L (5) 

The inequality Eq.(5) gives two regions e e z(q, 0) > 1 and e e ;(q, 0) < 0. The former case is commonly seen for most 
materials, such as metals where e(0,0) diverges positively, while the latter case corresponds to overscreening. Let us 
consider a positive test charge density ptest(o)- The test charge induces the screening charge density p SC reen(q) to 
reduce the energy of the system. The screening charge is always negative. Using this charge density, l/e e ;(q, 0) is 
expressed as 

' _ |Pscreen(q)| ^-g^ 



e ei (q,0) Ptest(q) 

When |/O SC reeri(q)| is larger than ptest(o) in Eq.(6), l/e e j(q, 0) is negative. This is the static overscreening effect. 

The microscopic origin of the negative dielectric function may come from the Following mechanism. The electric 
polarizability in covalent solids is qualitatively different from that in simple ionic crystals, since the covalency con- 
tributes to charge transfer between ions fl33| . For instance in ferroelectric oxides the polarization due to charge transfer 
is as large as the ionic polarization. In BaTiC>3 the nominal valence of Ti is +4 and the nominal d state configuration 
is dr. But the d-orbital of Ti and the p-orbital of O are strongly hybridized. This charge transfer produces current, 
and thus electronic polarization, which adds to the ionic polarization, making the effective valence (Born effective 
charge) of Ti twice as large |54],[55| ■ In undoped cuprate the same transfer occurs from the filled p- level of O to the 
empty upper Hubbard band of Cu, contributory to an extra polarizability. 

However, in a doped cuprate the situation is drastically different. The doped holes occupy mostly the oxygen 
p-levels, and they move either to the lower Hubbard band or the filled d z 2 orbital of Cu. Thus holes are transferred 
from O to Cu, creating the larger polarization, which adds to two ionic polarization. Then, Eqs. (1,2) give a negative 
dielectric function if x e i(q, 0) > 1. The more detailed discussions of the lattice softening using Born effective charge 



will be given elsewhere 36 



2o(q) =^To(q) + — rrz — * ?Tv\ — ' ( 7 ) 



While the in-plane Cu-0 bond-stretching LO phonon mode in the cuprates is strongly softened by doping, other 
phonon modes are relatively insensitive to the doping level |OJlj|jl6| . Thus, for the sake of simplicity, we assume that 
the basic lattice dynamics is not affected by doping, and the frequency of the bond-stretching mode is renormalized 
only by the interaction with charge fluctuations. This assumption allows to write the frequency corresponding to the 
maximum value of the spectral intensity of the renormalized LO phonon in terms of the electronic dielectric function 
as @|l|, 

^2 ,_ A 2 , ^£ o(q) 2 - ^Tp(q N2 
4(q> w 2o(q)) 

where luto and lolo are respectively the bare TO and LO phonon frequencies in the insulating state, and £gj(q, uj* lo ) 
is the real part of the electronic dielectric function. The experimental results that the frequency w£ (q) 2 of the LO 
mode is lower than that of the TO mode. It indicates that the real part of e e i (q, w£ (q)) is negative. In YBa2Cu306.95 
(YBCO) and Lai.85Sro.isCu04 (LSCO) the inversion of the LO/TO frequencies is most pronounced in the region q x 
= 0.25 to 0.75 and q y — -0.2 to 0.2, in the units of 27r/a, at oj* lo being the measured LO frequency (approximately 55 
meV for YBCO and 70 meV for LSCO) |l^,|l^,^|. We interpret this phenomenon as the consequence of the negative 

<./(q- qii- 

The quasi-particles are renormalized due to the on-site Coulomb interaction as mentioned in Introduction. The 
interaction increases the effective mass of the quasi-particles, but does not change the charge e of the particles. 
Therefore, the effective potential acting between quasi-particles k and k' is written as 

v «>^=^B)> (8) 

where V(q) is the bare Coulomb interaction and q is k— k'. Therefore with the normal screening the effective potential 
is always smaller than the bare potential. However, in the case of overscreening e' el (q,u>) can be negative, and thus 
e'<(q, w) contributes to the effective interaction and works to enhance the phonon mediated attractive interaction 
as seen in section III. Since the overscreening effect comes from various kinds of the correlation effect, e(q, ui) is a 
complicate function of charge, spin, and lattice. This effect is the core of the present mechanism. 

Consequently the attractive interaction should exist between the quasi-particles with k and k' when q = k — k' is 
in the q regions where the phonon softening occurs, as seen in Eq.(7). In addition the negative electronic dielectric 
function turns the repulsive electron-electron Coulomb interaction into attraction, as in the so-called negative-U 
mechanism, and can promote pairing. 
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III. FORMULATION FOR SUPERCONDUCTING PAIRING 



In this section we derive the equation to calculate the symmetry of the superconducting order parameter and T Cl 
using the effective interaction Eq.(8) 1 30 3]j| . In the effective interaction, the dynamical dielectric function e(q, uj) is 
given by the sum of the electronic dielectric function e e z(q, uj) and the ionic dielectric function e,; orl (q, uj) as 

e(q,w) = e e; (q, uj) + e ion {q, uj) - 1 (9) 

A minus unity in the right hand side of Eq. (9) comes from the fact that all the dielectric functions should be unity at 
the high frequency limit. We express the ionic dielectric function on a conventional form, 

ei „(q,a,) = ^-i° , (10) 
UJ Z — uj^o 

lulo and uito being respectively the frequencies of the bare longitudinal and transverse optical phonons in the 
insulating state. For simplicity, we consider one optical phonon mode, which seems to be most relevant to the 
superconductivity, and assume that it is dispersionlcss in the unrenormalized state. The spectral intensity function 
of total charge fluctuations is expressed by using Eqs.(4),(9-10) as 

p(q><*>) = -~Im[- 



7r e(q, uj) 



ij m[ 1 i + ^o(q)-^- 2 

* £ei(q^) e el (q,uj* LO ) 



where the first term in the right hand side of Eq.(ll) is the electronic spectral intensity p e i(q,oj), and u)* LO (q) is the 
LO phonon frequency renormalized by charge fluctuations and is given by Eq.(l). If we use the spectral representation 
for l/e(q, uj), the effective interaction Eq.(8) is written as 

V eff (q,oj)=V(q)/e(q,oj) = V(q)[l-2 frfn " p(q,n) ], (12) 

Jo Si — (U! + id) 

Using Eq.(12) we set up the Eliashberg equation linearized with respect to the gap function A(k, ito) as 

A(k, ioj n ) = -T J2 E V °ff ( k ~ k '> ^" ~ ™t) |r3^T . (13) 
i v w + ^ 

where u> = (2n + 1)ttT with n being integer, and is the quasi-particle energy measured from the Fermi level. We 
use an approximation that the damping of the quasi-particles is neglected. However, the modification of the band 
structure due to the correlation effect is taken into account by using the model band structure determined by the 
experimental results of angle-resolved photoemmision. If wc introduce the pair function defined by 

F(k,«j n ) = A(k,iw n )/(«2 + ££), (14) 

Eq.(13) is rewritten as 

F(k, iuj n ) = -—±-—tY^ E ^//( k - k '' lLU " ~ «^)F(k', ioj e ). (15) 
We then introduce again a function f{k,v) |37| defined by 



F(k,tuj n ) = 2 f 
Jo 



'w ^/( k ^) 

dv^r- (16) 



and also 



$(k) = 2|&| / dvf(k, 
Jo 



(17) 



Then, from Eqs.(14) and (16) we can show that $(k) is equal to i?eA(k, to a good approximation. The equation 
for <E>(k) is obtained substituting Eqs.(16) and (17) for Eq.(15) as 
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*(k) = -^^(k,kO tan ^ k ' /2T ^(kO, (18) 

In Eq.(18), the kernel K(k, k') is given by 

f°° n(U - k' CI) 

If we insert Eq.(ll) for p(q,ui) in Eq.(19) and use the Kramers-Kronig relation for 7m(l/e e /(q, ui)) in p e ;(q,w), the 
kernel Eq.(19) is written as 



TO 



K(k 7 k')-v(k k')[ £e;(qj0) +2^ dn m + ] ^ + M p el (c l ,n) £ei(qj ^ o(q))w2o(q)Ko(q) + | a | + | 6 ,| r 

(20) 

Roughly speaking, in the right hand side of Eq.(20) the first two terms are the electronic contribution and the third 
term is the phonon contribution. However, the electronic contribution is mixed also in the third term as seen in 
Eq.(20). Therefore, the kernel has a vibronic nature. 

The spectral intensity p(q,u) is approximately written around u> being w* LO {q) ]30[ | as 

n(a . o _ 1 ^lo- "to A (q) r , n 

m ' ' <,(q,u£ (q))» 2™j (q) [ W - ^(q)]* + A*(q) 1 J 

Afn x = ^(q^Lo) ^Lo(q)-^To(q) f22 ^ 
w <i(q^£o(q)) 2 2c£ ( q ) 1 j 

In Eq.(22), e" ; (q, w£ ) is the imaginary part of the electronic dielectric function. A(q) is the half width in the 
frequency u* LO °f ^ ne spectral function. If we use the experimental frequency values in Fig. 2, e' e; (q, lo* lo ) is given 
-0.068 from Eq. (7) for YBCO. Using this value, Eq. (22) and the neutron scattering experimental value of A(q) ~ 
2meV, we obtain e"j(q, ijJ^o) *° ^ e 0-0069. Therefore, we can approximate e e ;(q, w£ (q)) by e' eZ (q, Cc>2 ), while we can 
see that the damping on the low frequency charge fluctuation is very weak. 

In conventional systems with ^(q, w) > 1 the third term in Eq.(20) cannot be enhanced. However, when ^(q, w) < 
and its absolute value is small, the phonon contribution can be strongly enhanced, as seen from Eq.(20). The strong 
softening of u>* LO (q) also enhances the third term. 



with 



IV. QUASI-PARTICLE BAND STRUCTURE AND ELECTRONIC DIELECTRIC FUNCTION IN 

CUPRATE SUPERCONDUCTORS 

In this section, the quasi-particle band structure and the dielectric function are determined by utilizing the ex- 
perimental data of the angle-resolved photoemission and inelastic neutron scattering. The results will be used in 
the following section in calculating the symmetry of the superconducting order parameter and the superconducting 
transition temperature with the Eliashberg equation. 



A. Quasi-particle Band Structure 



We assume that the layers responsible for the high T c superconductivity are mainly the CuC>2 layers. Then, we 
construct the energy band of the CuC>2 layer (the a — b plane) as to agree with the results of the angle-resolved pho- 
toemission spectroscopy (ARPES) |2£]j&| . The band structure is almost universal for all the cuprate superconductors 
with optimum doping |2f| . The angle-resolved photoemission j}8) and the calculation including the correlation effect 
HI shows the energy band in the cuprate superconductors is very flat in the region of k x (k y ) = 0.25 — 0.75 in units of 
27r/a, at 10 — 30 meV below the Fermi level. Hereafter, we express the wave number in units of 2ir/a. This is the region 
where the superconducting gap is the largest, and the kink in dispersion, which is interpreted as an evidence of strong 
electron-phonon coupling is most conspicuous. Therefore it is reasonable to expect that the charge fluctuation 
created by quasiparticles on the flat band is very low in energy, comparable to the optical phonon frequency. Then, 
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the quasi-particles and phonons can be highly mixed, and this system consequently becomes vibronic. As a result the 
quasiparticle-phonon interaction is highly enhanced, as will be shown in sections III and V. 

The band structure with the flat regions along the k x and k y axes which originate from the strong renormalization 
effect ]3i| is reproduced by the following function, 

Cfc = Cs(to + 2ti(cos2Trk x + cos2nk y ) + At2Cos2Ttk x cos2Ttk y + 2t^{cosA j nk x + cosAirk y )) + Fg(Ef — uj g ) — Ef, (23) 

where to, ti, £2, and t% are respectively 0.0, —0.2, 0.0 and —0.04, and Ef is —0.11 and oj g is from 0.01 to 0.03 in units 
of eV. In this equation,if the second term is dropped, the equation approximately coincides with that of the usual 
LDA band calculation Q . The second term is required to reproduce the flat regions mentioned above. The function 
Fb is given by 

F B =W X (1+ ^tanh((k x - k L )/X L ) - ^tanh({k x + k L )/\ L )) + W y {l + ^tanh{{k y - k L )/\ L ) - ^tanh{(k y + k L )/X L )), 

(24) 

with 

W x = -^tanh{(k v - k T )/\ T ) + ^tanh((k y + hr)/X T ), (25) 

W y = -^tanh((k y - k T )/\ T ) + ^tanh((k y + k T )/\ T ), (26) 

where k^, kx and Aj, = are the parameters to control the extent of the flat regions. A function Fb has sizable 
values only in certain regions near the k x and k y axes. For Cb and Fb, a constraint Cb + Fb = 1 is imposed. Fig. 
1 shows the band structure calculated by using Eqs. (23-26) where fcj,, kr, and Xl(= At) are, respectively, taken as 
0.16, 0.15, and 0.2 in units of 2ir/a. The values above are chosen to reproduce the ARPES data for the optimally 
doped YBCO. 

As seen in the inset (a) of Fig. 1 the flat regions of the energy bands are located along the k x axis and extends 
from (-0.5, 0) to (-0.25, 0) and from (0.25, 0) to (0.50, 0), and similarly along the k y axis. The behavior has been 
observed in the optimally doped YBCO and LSCO by ARPES J25|,|3^]. The energy of the flat band is from 10 meV 
to 30 meV below the Fermi level (^5|. In Fig. 1, uj g is chosen to be 18 meV below the Fermi level. These flat regions 
mainly originate from the correlation effect reflecting the electron scattering due to spin fluctuations [^9[ . In these 
regions the effective masses and the damping constants of the quasi-particles are very large. On the other hand, the 
energy band in the directions from the (0, 0) point to the (0.5, 0.5), (0.5, -0.5), (-0.5, 0.5), (-0.5, -0.5) points is not 
much affected by the correlation effects and the dispersion of the electronic energy band structure is usual as seen in 
the inset (b) of Fig. 1. 
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(a) (b) 




FIG. 1. An typical example of electronic band structure for calculating Tc and the superconducting order parameter. The 
parameters required for giving the flat band region as seen in this figure are fez, = 0.16, kr = 0.15, and Al(At) = 0.2 in unit 
of 2-w/a, and the parameter for the depth of the flat region uj g — 0.018eV\ The z-axis indicates the energy of quasi-particles 
in units of leV. (a)The cut of the band structure from (—0.5,0) to (0.5,0). The Fermi level is displayed, (b) The cut from 
(-0.5,-0.5) to (0.5,0.5). 



B. Dielectric Function 



Recent inelastic neutron scattering measurements on the optimally doped YBCO show that the frequencies of the 
LO and TO phonons of Cu — O bond-stretching LO mode are inverted in the region of q x = 0.25 to 0.5 [16 19|,|20|,p6]| , 



suggesting that e! el (q x ,uj*) in this region is negative according to Eq.(7) and the LO phonons are overscreened by 
charge fluctuations. It was also observed that the phonon dispersions are anisotropic along the a— and 6— axes. In 



7 



this paper, however, we neglect this anisotropy for the sake of simplicity, and assume the tetragonal symmetry of the 
crystal skipping details of the dispersion. The effect of the anisotropy will be discussed elsewhere jyj . 



CO 



CO L0 7 3meV 
CO T0 7 2meV 



CO L0 55meV 



q 



.25 



0.75 



X 



(q y ) 



( X27i/a ) 



FIG. 2. The schematic view of the difference in the energy levels of the breathing LO and TO phonons between un-doped 
and doped cuprates for YBCO. In the undoped state, the LO and TO branches are dispersionless, and are at 73 and 72 meV, 
respectively, while in the doped state the LO level goes down to 55meV in the range from 0.25 to 0.75 in q x and q y directions 
although the TO level still remains dispersionless. The q x (q y ) is measured in units of n /a 



We then assume that the dielectric anomalies occur in the shaded regions as shown in Fig. 3. We divide the q 
space into two regions, Region II is the shaded area in Fig. 3, where anomalous dielectric behavior is observed, and 
Region I which is outside Region II with normal dielectric behavior. In Region I we assume the Thomas-Fermi type 
of screenings of the Coulomb interaction due to quasi-particles, since the plasma frequency is much higher than kT c 
and the frequency of the optical phonons. Then, the inverse of the dielectric function in Region I is expressed as 



1 



1 



e ei (q,w) e e ;(q,0) 



2 : 
QTF 



(27) 



where qxF is the Thomas- Fermi wave number and its value is an order of unity in the units of the reciprocal lattice 
constant for the optimally doped cuprate superconductors JM2| . In this region, we note that the repulsive interaction 
by the spin fluctuation scattering is located at (0.5,0.5) f|,|7|,|8| , but in the present paper this interaction is neglected 
in order to concentrate on the dielectric charge channel. The repulsive interaction actually exists, but for the moment 
we presume that the main pairing contribution arises from the dielectric charge channel. In Region II, we assume the 
excitation frequency of the electronic charge fluctuations is much lower than the plasma frequency for q in the region. 
The existence of such a low energy charge fluctuation due to anomalous electronic structures has been suggested in 
P3| . The existence of the low frequency charge excitation is a necessary condition for the appearance of the negative 
dielectric function. Then, in Region II, for the inverse of the dielectric function for positive low frequencies, we assume 
the following one pole approximation 
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_J_~A( -), (28) 

e e j(q, (jj) u-UQ + i-y 

where A, u>q, and 7 are respectively the amplitude, the frequency of the charge fluctuation density maximum, and the 
damping constant. From Eq.(28), A is given by 

A=^\ (29) 

Using Eq.(7) and the experimental values of the frequencies shown in Fig. 2 we estimate ^(q, <~o* LO ) to be -0.068. The 
value of ujo is much lower than the plasma frequency {%oj p = 0.8eU) measured at small q by reflectivity. Inserting 
the value of w£ for ui in Eq.(29), we obtain A to be from 132 meV if hui n is assumed to be 64mcV. The value of A 
linearly increases with luq. However, we note that Eq.(28) is valid only in the vicinity of the low charge excitation 
energy hui . Using the value of A and e" ; (q,w£ ) calculated from Eq.(22), we obtain 7 to be ~ 1 meV. 
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FIG. 3. The schematic view of the division into the region I and the region II (the hatched region) in the space of the 
scattering wave number between two quasi-particles. In the region I, the Coulomb interaction is screened as in a usual metal, 
while anomalous screening occurs in the region II due to strong correlation effects. Thus, for the dielectric function, the 
expression for Thomas-Fermi type of screening is used in the region I, while the one obtained from the experimental result 
(neutron scattering) is used for the region II. 



V. SUPERCONDUCTING GAP SYMMETRY AND T c 



In the previous section, we constructed the model electronic band structure and derived the dielectric function 
based on the recent experimental results. In this section, we determine the kernel, Eq.(20), using the experimental 
values given in section VI and calculate T c and the superconducting gap function by solving Eq.(18). 
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Let us first concentrate on the kernel function in Eq.(20). The first term in the bracket in Eq.(20) represents the 
static dielectric function. It gives the repulsive contribution in the region I in Fig. 3, while the term becomes negative 
causing the attractive interaction in the region II. The term in the region II is determined by extrapolating the real 
part of the dielectric function Eq.(27) to the zero frequency as 

' -A— (30) 



4(q,0) uj 

where Eq.(28) is used. At this moment we do not have exact information on the value of u>q. However, since strong 
phonon-induced charge transfer is expected, the value of luq must be close to the value of lo* lo , 55 meV. We thus 
assume the value of ojq to be 64 meV. Thus, l/e' el (q, 0) is given to be ~ 2. The second term always vanishes on the 
Fermi surface since £fc = = on the surface, but it gives the repulsive interaction in finite energy ranges from the 
Fermi level. It is found from the numerical calculations that the contribution of the second term dominates over the 
attractive first and third term when and are sufficiently far from the Fermi level (4]J| . The third term gives the 
attractive interaction originated from the optical phonons which is renormalized by the electronic dielectric function. 
This term is strongly enhanced due to overscreening in the region II. Such a strong enhancement in the region II can 
understood by rewriting and comparing the third term for region I and II, respectively, as follows, 



and 



i 2 ^o 2 (q)-^ro 

q 2 + q^ F uj* L0 2 (q) 



1 ^* LO (q) - w* Q 

e/ a i( < J. w £o) <4o 2 (<i) 



for the region I, (31) 



for the region II, (32) 



where both and in the third term are set to be zero to compare their contributions on the Fermi surface. If q is 
set to be 0.5 for both equations, the reasonable value of qrF is employed and the observed values uj* lo in both regions 
are substituted, it is found that the contribution in the region II becomes more than about 100 times larger than 
that in the region I. The reason is explained as follows. In the region II, the amplitude of l/e*,(q, u>) almost diverges 
negatively at oj - Therefore, the amplitude of l/e' ei (q, to) at u> w w£o s ^ n remains very large as shown in Fig. 4. The 
value of Eq.(32) is estimated to be 10.7. On the contrary such a anomalous structure does not exist in the vicinity 
of the phonon frequency in the region I and the amplitude of l/e^(q, u>) remains normal values as in conventional 
metals. These different features for the region I and II are seen in Fig. 4. 
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FIG. 4. The uj dependences of l/e' el {oj) in the region I and II and — Im[l/e e i(u))] in the region II. The solid green line represents 
l/e' el (u>) in the region I, while the solid black and the dashed red lines stand for l/e' el (uj) and — Im[l/e e i(uj)], respectively. In 
the region I there is no characteristic peak structure within the frequency range (0 ~ O.lel^) because its peak is located at the 
high plasma frequency range (« 0.8eV). On the other hand, in the region II the fluctuation peak lies at ujq as shown in this 
figure. 
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Thus, the electric overscreening effect is found to enhance the phonon contribution strongly and to give a strong 
attractive interaction. In fact, the comparison between the first and the third term contribution in the region II for 
the quasi-particles on the Fermi surface tells us that the third term contribution is 10 times larger than the first term 
in this case. This result indicates that the electron-phonon interaction enhanced by the overscreening could be the 
main pairing interaction in the cuprate superconductors. 

Next, we actually calculate T c and the k dependence of the gap function by solving Eq. (18). For this purpose we 
use the following technique. We introduce an eigenvalue equation with an eig envalue \{T) 

A(T)*(k) = -^mX) tanh ^' /2T) HX). (33) 

We numerically calculate A(T) as a function of temperature T. The temperature satisfying A(T) = 1 corresponds 
to the superconducting transition temperature T c . In addition, since $(fc) corresponds to the superconducting gap 
function, the /e-dependence of $(fc) gives the superconducting gap symmetry. However, it should be noted that the 
amplitude of the gap function is meaningless since the equation is a linearized equation which is valid only in the 
close vicinity of T c . 




FIG. 5. A typical exmaple of the superconducting gap function in the wave number space. The x and y directions indicate k x 
and k y , i.e., quasi-particle wave number, respectively, and the z-axis gives the amplitude of the gap function. In the numerical 
calculation, the parameters ui g , uio, Kl, Kt, and Al are set to be 13meV, 64meV, 0.16, 0.15, and 0.2, respectively. The used 
quasi-particle band structure is displayed in Fig. 1. 



A typical example of the superconducting gap function is shown in Fig. 5, where we use 13meV for uj g in Eq.(23) 
and 64meV for wo, and the other control parameters are the same as those employed to depict the band structure 
shown in Fig. 1 which corresponds to the optimally doped case. The calculation in this case shows that T c is about 
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180K. As seen in Fig. 5 it is found that the gap function has c4 2 -fc 2 wave symmetry. The d k 2_ k 2 symmetry is possible 
under not only the spin fluctuation mediating pairing but also other anisotropic pairing interactions j44| . We note 
that the gap function has nodes while it shows finite amplitudes even at the center of the flat region in the electronic 
structure. This means that the attraction is so strong that the gap function does not vanish even in regions except 
for the Fermi surface. This fact is consistent with the experimental results in ARPES for the optimally doped case 
showing high-T c p5| . Here, let us explain why the d k i_ k i symmetry appears in this case. The gap amplitude is 
generally dependent on the number density of the interacting quasi-particles and the strength of attraction. Since the 
number density of the quasi-particles in the flat region is much larger than other regions, the gap function is expected 
to grow considerably in the flat region. A pair of quasi-particles on the opposite sides of the same flat regions along 

tliG kx (fey ) 

axis are attracted, while the interaction is repulsive for a pair of quasi-particles belonging to the flat regions 
that are rotated by 90° relative to each other. Thus, it is found that the gap amplitudes on the different axes should 
show different signs, resulting in the d-symmetry. 
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FIG. 6. The w s (eV) dependence of T c (eV). The other paramters are the same as ones used in obtaining Fig. 5. 



As seen in Fig. 6, we found that the value of T c is nearly inversely related to the depth of the flat region. This 
tendency is consistent with the experimental results of ARPES j2j| , that the depth of the flat region increases with 
increasing doping, and crosses the Fermi level in the overdoped range of composition. At the optimum doping, the 
flat region shows a minimum depth. On the other hand, the value of T c is dependent also on the width Wu of Region 
II, increasing almost linearly with the width. Our calculations show that T c exceeds 300K when the value of Wu 
approaches to 0.4 with fixing tv g to be 13meV, corresponding to the optimum doping. It is quite encouraging that 
such a high value of T c can be achieved with the present model. 



VI. DISCUSSIONS AND CONCLUSION 

In this paper the mechanism of high T c superconductivity based upon overscreening of phonons |3(],|3l| was developed 
further using the results of two recent experimental observations. One is the neutron scattering measurements showing 
that the frequencies of a bond stretching LO phonons in hole-doped superconducting cuprate oxides being strongly 
softened in the wave number regions around the Brillouin zone boundaries along the k x and hy and that the frequency 
of the LO phonon is much lower than that of the TO phonons in these regions [|l6| , ^9|j20[^r|] . This phenomenon 
is explained in terms of the overscreening of the ionic polarization associated with the LO phonon by phonon- 
induccd charge transfer that results in the negative electronic dielectric function. Another one is the angle-resolved 
photocmission measurement indicating that in the optimally hole-doped cuprates the flat energy band appears just 
below the Fermi level around the X (M) points in the Brillouin zone |25| , |3"8| . Using the pair interaction derived from 
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the neutron scattering results and the band structure mentioned above, we set up the Eliashberg equation pO| , pi| . 
Solving the equation, we observed that the symmetry of the superconducting order parameter is of d^-k 2 and the 
transition temperature is in excess of 200K. 

The most dominant contribution to the superconducting pairing in the kernel of the Eliashberg equation originates 
from the phonon contribution that is strongly enhanced by the overscreening effect of the Coulomb interaction. Up 
to now, several phonon mechanisms leading to d-wave superconducting gap function due to the anisotropy of the 
electronic structure have been suggested J45j. In contrast, the present theory relies upon a mechanism of phonon 
overscreening resulting in the negative electronic dielectric function and strong enhancement of its attractive contri- 
bution. As the most of the phonon mechanisms, the present theory predicts that the isotope effect appears in 
the cuprate superconductors with a relatively low T c such as LSCO, but this effect may not be visible in the cuprates 
with a high T c such as YBCO and Bi-2212, as observed by experiments. This result comes from the fact that as the 
superconducting transition temperature increases, the pure phonon contribution for the pairing diminishes and the 
electronic contribution increases instead. This tendency is consistent with the experimental results |Q. 

The softening of the LO phonon near the Brillouin zone boundary was observed by neutron scattering even in 
the superconducting Bao.6Ko.4Bi03 with negligible spin fluctuation |47j], and non-superconducting Lai.69Sro.3iNi04 
pif and Lao.7Sro.3Mn03 |49|. It appears that the overscreening effect is a common nature in hole-doped transition 
metal oxide systems. This phenomenon may be explained in terms of the "negative Born effective charge", which will 
be further elaborated later This concept nicely explains the phonon softening observed for the Peiels-Hubbard 
Hamiltonian in one dimensional Cu-0 chain model [pr| , ^8| . 

In the cuprates the bare bandwidth, calculated for instance with the LDA, is of the order of 1 eV. Band-narrowing is 
brought about by spins that almost localize charges and create the extended saddle point in the electronic band struc- 
ture. In the present paper this effect was considered phenomenologically by using the renormalized band structure. 
However, a more satisfactory theory should include the quantum effect of spins more explicitly. Using the p — d model 
hole doping was found to create in-gap states in the Hubbard gap just at the Fermi level, with a dominant oxygen 
p-character J50|. A better model band structure with these features needs to be developed. A problem here is that 
we have much less information about the unoccupied states compared to the occupied states, since the photoemission 
experiments provide information only for the latter, while the inverse-photoemission technique, which should give the 
information on the former, has much less resolution. 

As we mentioned above in the present mechanism LO phonons excite holes from the local oxygen p-state to the 
Cu d-state. The results of inelastic neutron scattering measurements on YBa2Cu3 0g.95 [pl[ suggest the final state of 
this transition could be the d z 2 level of Cu, since strong mixing of the in-plane Cu-0 bond-stretching mode and the 
apical oxygen mode (62 meV) was observed. This makes the two-band phononic model [(52) relevant in the present 
context. In addition there is a possibility of the transition from an oxygen in the rr-direction to another oxygen in the 
y-direction contributing to the e — p coupling (53). In the present calculation just one band was assumed for the sake 
of simplicity, but the reality is likely to be more complex. 

Another point that has not been included in the present paper, and has to be addressed in future publications, is the 
question of the in-plane anisotropy and the stripe fluctuations. Both the inelastic neutron scattering measurement of 
the phonon dispersion and the photoemission measurement detected surprisingly strong anisotropy in the Cu02 plane 
p9| , p0 54 3q ]. This anisotropy most likely is related to the spin-ch arg e phase separation in the form of stripes [ p5[ . 
While static stripes apparently compete against superconductivity p7| , it is reasonable to assume that the propensity 
for stripe formation remains strong even in the superconducting phase |5^,^,^,^|. The presence of stripes or 
stripe fluctuations will produce further narrowing of the band and spin-charge interference phenomena, which should 
enhance the superconductivity. The possibility of such enhancement by confinement, in a more general sense, was 
emphasized by Phillips |58[ . However, since we do not have sufficient information to formulate these phenomena into 
a rigorous theory, in the present paper we have not explicitly taken them into account. Further experimental as well 
as theoretical researches are needed to fully account for the spin-charge synergy. 

In the present mechanism the role of spins is merely to bring down the energy scale of holes to the level of phonons, so 
that the vibronic resonance can take place. We do not, however, exclude the possibility that the magnetic mechanism 
also contributes to pairing, at least in some ranges of composition. Even in such a case the phonon mechanism does not 
compete against the magnetic mechanism because of the d-symmetry. It is most likely that the relative importance of 
the phononic and magnetic mechanisms depends upon charge density; magnetic mechanism could be more important 
in the underdoped region, while the phononic mechanism may dominate the optimum and overdoped region. 

In the present paper we sacrificed some details as discussed above and greatly simplified the model. Furthermore it 
is well known that the calculation of T c in the Eliashberg theory strongly depends upon details of the parameters. The 
purpose of this work, therefore, is not to try to describe the superconductivity in the cuprates very accurately, but 
to demonstrate the possibility of the phononic mechanism of high-T c superconductivity based upon the overscreening 
phenomenon. The estimated value of T c is even higher than the experimental values over wide ranges of parameters. 
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This is in part due to the fact that the calculation of Tc is essentially based on a mean field theory and does not 
include the superconducting fluctuation effect. Inclusion of the fluctuations should significantly decrease T c , since the 
pairing force is strong and the system has a two dimensional nature. The incoherent superconducting fluctuations 
will be observed as the pseudo gap above T c . 

While there are a number of limitations in the present model because of the simplifications and the phenomenological 
nature, the result described here strongly suggests that the phononic mechanism has to be taken seriously, and careful 
studies are warranted. We will address some of the problems that are neglected here in future publications. 
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